Astronomy Sz Astrophysics manuscript no. Contursi2007 


© ESO 2008 


February 5, 2008 





ISOCAM Mid-infrared spectroscopy and NIR photometry of the 

HI I complex N4 in LMC. * 

A. Contursi^ M Rubio^, M. Sauvage^, D. Cesarsky^'^ R.Barba^, and F. Boulanger^ 



1 



Max-Planck-Institut fiir Extraterrestrische Physik (MPE), Postfach 1312, 85741 Garching, Germany 
^ Departamento de Astronomi'a Universidad de Chile, Casilla 36-D, Snatiago, Chile 

^ CEA/DSM/DAPNIA/Service d'Astrophysique, CE Saclay, 91191 Gif-sur-Yvette Cedex, France, France 
^ Departamento de Fi'sica, Universidad de La Serena, Chile 

^ Institut d'Astrophysique Spatiale, Bat. 121, Universite Paris XI, 91450 Orsay CEDEX, France 
Received 22 August 2006/ Accepted 

ABSTRACT 

Context. We present the analysis of ISOCAM-CVF and J,H and Ks photometry data of the HIT region complex N4 in 
the Large Magellanic Cloud (LMC). 

Aims. The aim is twofold: 1) to study the connection between the interstellar medium and the star content of this 
region; 2) to investigate the effects of the lower than galactic metallicity on dust properties. 

Methods. A dust features - gas lines - continuum fitting technique on the whole ISOCAM-CVF data cube, allows the 
production of images in each single emission and the detailed analysis of dust (both continuum and bands) , and ionized 
gas. The Near Infrared (NIR) photometry provides, for the first time, information on the stellar content of N4. 
Results. The Mid-Infrared (MIR) spectral characteristics of N4 are those expected for an HII region complex, i.e. very 
similar to those observed in galactic HII regions. The images in single dust feature bands and gas lines clearly show 
that the HII region core is completely devoid of the carriers responsible for the Aromatic Features (AFs). On the other 
hand, the ionized gas arises almost completely in this dust cavity, where also the two main exciting stars of N4 are 
located. HII region models from Stasinska (1982) predict an HII region size which corresponds to the observed size of 
the dust cavity. 

We find evidences that the effect of lower than Galactic metallicity (although not extreme as in the case of LMC) on 
the carriers responsible for the AFs, is not to prevent their formation or to modify their chemical properties, but to 
enhance their destruction by the high and hard interstellar radiation field. We argue that this is the dominant process 
responsible for the absence of AFs in the HII region core. We show that this mechanism is more efficient on smaller 
dust particles/molecules thus affecting the dust-size distribution. We argue that effects on dust-size distribution, rather 
than the different dust properties due to a lower metallicity, should be taken into account when analyzing more distant 
relatively low metallicity galaxies. 

Finally, the analysis of the stellar content of N4 reveals 7 stars: 4 reddened O MS stars and three stars with envelopes. 
In particular, one of these, seems to be an Ultra Compact HII region containing an embedded Young Stellar Object. 

Key words, interstellar medium: HII regions - dust features - gas lines - continuum 



1. Introduction 

N4 (Henize I1956P is an HII complex in the north-west 
part of the Large Magellanic Cloud (LMC). The physical 
characteristics of the ionized gas and associated molecular 
cloud, have been extensively studied by Heydari-Malayeri 
and Lecavelier des Etangs (|1994p through optical and sub- 
millimeter line spectroscopy. 

The Ha emission of N4 (Figure [1] kindly provided by 
Heydari-Malayeri) clearly shows that this HII complex is 
composed of two ionized nebulae: N4A and N4B. In this 
paper we will concentrate only on N4A. This is brighter 
and younger than N4B. It is composed of a bright dense 
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front to its north-east part and a more diffuse low sur- 
face brightness, partially density bounded, component to 
its south part. The main ionizing sources in N4A are two 
stars not visible at optical wavelengths (marked in Figure 
[l]as A and B). Heydari-Malayeri and Lecavelier des Etangs 
(|1994p calculated that the main ionizing sources correspond 
to either one 6OM0 star or two 40 M0 stars. 
The average extinction in N4 derived from optical gas emis- 
sion is quite low, Av'=0.4 although it is presumably higher 
in the north-east front where the gas density reaches its 
maximum. The parental molecular cloud is formed by two 
spatial components: one peaking on N4A and the other 
peaking to the east of the HII region. The molecular cloud 
has two velocity components: a high velocity optically thick 
red component and a thin blue molecular sheet. This sheet 
is probably pushed towards us by the HII region pressure 
(Heydari-Malayeri and Lecavelier des Etangs fl994p . 

A complementary study on the hot dust emission in 
N4A has been done by Contursi et al. (|1998p using broad 
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band ISOCAM filters centered at 7 and 15 /im. At these 
wavelengths, the emission is dominated by the Aromatic 
Features seen in Emission (AFEs) between 3 and 10 /im and 
a continuum steeply rising at wavelengths longer than ~10 
lim. The overall Mid-Infrared (MIR) morphology in N4A 
is generally quite similar to the CO emission, although the 
MIR and CO peaks do not exactly coincide. The strongest 
MIR emission comes from the dense north-east front as ex- 
pected, since here the density of the matter is the highest. 
The 15/7 /xm ratio peaks in the region of N4A where the two 
main exciting stars are supposed to be. Its high value, typ- 
ical for an HII region, indicates that the emission at 15 fim 
is dominated by a strong continuum and that the emission 
at 7 fim is significantly reduced. Contursi et al. (|1998p 
suggested that this decrease is primarily due to destruc- 
tion of the AFEs carriers in hard and strong Interstellar 
Radiation Fields (ISRF). 

This analysis was part of a wide ISOCAM campaign 
aimed at studying the interstellar medium (ISM) prop- 
erties in the Magellanic Clouds with particular emphasis 
on the role played by metallicity, taking full advantage of 
the unprecedent (at that time) spatial resolution provided 
by ISO at these wavelengths and to the proximity of the 
Magellanic Clouds. One of the outcomes of the broad band 
ISOCAM imaging study of N4 was that this region appears 
as a classical Galactic HII complex, composed mainly of 
three phases: 1) The HII region itself close to the exciting 
stars corresponding to the maximum of the 15/7/im ratio; 
2) the interface between the HII region and the molecu- 
lar cloud (Photon-Dominated Region or PDR) where the 
MIR emission reaches its maximum due to both high dust 
density and relatively high ISRF; 3) the parental molecular 
cloud, where the MIR emission is still dominated by AFEs 
probably arising from the outskirts of the molecular clouds 
and excited mainly by the general LMC ISRF. 

In this paper we present a more detailed study of the 
MIR emission properties of N4A analyzing a very high 
signal to noise spectro-imaging data set carried out with 
the Circular Variable Filter (CVF) on ISOCAM. Thanks 
to the spectral capabilities, we are able to directly show 
that our previous physical interpretation of the broad band 
MIR photometry data of N4, were correct. We also present 
new high resolution J, H and Ks images, taken at Las 
Campanas Observatory (Chile). All these data allow to bet- 
ter investigate the existence of detailed metallicity signa- 
tures on the dust emission which do not show up in the 
broad band emission and to relate together dust, ionized 
gas and stellar content. These data provide high resolution 
analogs for the study of the ISM in conditions similar to 
those of young galaxies in the distant Universe which are 
becoming more and more available with the advent of large 
ground based telescopes (VLT, Keck, ALMA) and space 
missions (SPITZER, HERSCHEL). 

2. The data 

2.1. ISOCAM-CVF data redution 

A complete scan of the two long-wavelength Circular 
Variable Filters (CVF) was performed going down in wave- 
lengths, step by step, first from 16.2 to 9.0/ito (LW-CVF2) 
and second from 9.4 to 5.0 fim (LW-CVFl). At each CVF 
step, between 7 and 11 frames of 2.1 sec were taken. The 
total observing time was 2974 sec. The pixel field of view 
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Fig. 1. Ha emission of the LMC HII complex N4A and N4B 
(from Heydari-Malayeri and Lecavelier des Etangs 1994). 
The image has been rotated to the same orientation as the 
ISOCAM-CVF images. 



was 6" and the full field of the observation 3' x 3' cen- 
tered on RA=4:52:8.2 and DEC= -66:55:16 (J2000). The 
data reduction was performed as described in Boulanger et 
al. (|2004p . The reduced d ata cube can be retrieved from 
the ISO archive at the site |http://www. iso.vilspa.esa.es as 
Highly Porcessed Data Products (HPDP) under the name 
Mid-IR Spectro Imaging ISOCAM CVF Observations. 

2.2. Near Infrared Broad band filters data 

Near Infrared J, H and Kg deep images of N4A were ob- 
tained with the Dupont 2.5m telescope at Las Campanas 
Observatory (Chile) in the night of December 27, 1996 with 
the 256 X 256 NICMOS III camera IRCAM (Persson et al. 
1992). The spatial resolution of the system was 0.35" /pixel, 
and the typical seeing was 0.9 "in Kg. 

The observation consisted of a series of 10 frames, each 
individual frame with an integration time of 20 sec per fil- 
ter. This procedure was repeated in a 9 position mosaic 
with separation of 20 ". The total integration time in each 
filter resulted in 1800 sec in Ks, 1800 sec in H, and 2000 
sec in J. Sky frames in each filter were taken in a field with 
faint stars and no extended emission at 12"W and 134"S 
of the source position. The sky field was observed in a sim- 
ilar way as N4. To produce the final images, each image 
was flat fielded and sky subtracted, and then median av- 
eraged combined using IRAF procedures. The final images 
were registered with respect to the Kg image by means of 
several common stars. The resulting mosaic covers a 65" x 
76"area. 

Aperture photometry in each filter was performed with 
IRAF/DAOPHOT on the mosaic. For the purpose of this 
study we concentrated on the sources brighter than Ks <13. 
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Fig. 2. ISOCAM CVF spectra of four different regions of N4A. The regions where the spectra were extracted are shown 
in each panel on the gray scale image in the 7.7 /im aromatic feature. Note that the steep rise or decrease at the very 
end of the spectrum are instrument artifacts and should not be considered as real. 



The stars associated to N4A are shown in Figure 5. Table 
1 gives the NIR magnitudes and colors of these stars. The 
photometric errors are axs =0.06 mag, an = 0.05 and 
a J =0.03. Photometry of stars in our mosaic far from the 
N4A nebulosity were compared to their 2MASS photometry 
showing a good agreement. We estimated the flux contri- 
bution from the nebulosity to our measured photometry to 
be less than 10% and 5% in the in the Kg J bands, respec- 
tively. 

3. Results 

3.1. The Mid-Infrared Spectrum of N4A 

Figure [2] shows the MIR spectra of 4 regions in N4A repre- 
sentative of 4 typical components of an HII complex. The 
upper left panel spectra (a) were produced by averaging 
4 pixels (6"/pix) centered in the position of the two main 
exciting stars in N4A. The upper right panel (b) shows 
the spectrum of the ionized gas peak emission. The bottom 
left spectrum (c) is the average of 18 pixels on the bright 
north-east front of N4A; the bottom right spectrum (d) cor- 
responds to the CO eastern peak, outside the ionized cloud 
(Contursi et al.l998). 

The spectra around the stars (a) and on the north-east 
front (c), are dominated by the AFEs and a continuum 



longward of 10 /im which becomes more and more promi- 
nent approaching the stars (b and a) . AFEs are stronger in 
the north-east front of N4A, as expected in a region where 
both density and ISRF are high. This spectrum (c) is typ- 
ical of PDRs confirming that this region is the interface 
between the HII region and the molecular cloud of N4A. 
In the spectra (panel a) corresponding to the regions close 
to the stars the prominent features are the fine structure 
gas lines of SIV at 10.5 /im, Nell at 12.8 fj,m (although 
blended with the AFE at 12.7 /im) and Nelll at 15.6 /im. 
The AFEs over continuum ratio reaches its minimum in the 
region where the fine structure line are the most intense (b) 
and AFEs practically disappear. The spectrum (d) towards 
the eastern CO peak presents very weak AFEs and almost 
no continuum. These are characteristics of quiescent regions 
and suggest that the emission arises from a thin external 
shell of the molecular cloud probably corresponding to its 
red velocity component. 

From this global spectral analysis it appears that N4 
MIR dust and gas properties are very similar to what is 
found in galactic HII complexes (Cesarsky et al. I1996b[ 
1995a|) . i.e. its MIR emission spectrum is typical of a 
HII/PDR/Molecular Cloud combination. 
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3.2. Emission in the single dust bands and gas lines 

Thanks to the very high signal to noise and the fact that 
we have spatial and spectral information on a region of 
3'x3' centered on N4A, this data set is perfectly suitable 
for a detailed spectral and spatial analysis in the single 
AFEs, gas lines and pure continuum. 

We produced images in each dust feature (namely, 6.2 //m, 
7.7 /im 8.6 /im 11.3 /im) and fine structure lines (SIV, 
Nell and Nelll) in the following way. Following Boulanger 
et al. 119981 we performed Lorentzian fit on the dust bands 
and Gaussian fit on the gas lines, to the whole data cube, 
producing maps in each dust feature and fine structure line 
of the ionized gas. First we fit together the 6.2 ^m, 7.7 
and 8.6 /Ltm bands and a straight line as continuum letting 
the slope vary at each pixel. Then, we fit together the 
SIV, 11.3 NeII(+12.7 ^m) and Nelll lines, and a straight 
line as continuum. The reason why we separately fit the 
two parts of the spectrum is because in the HII region 
core the continuum shortward and longward 10 /zm might 
not have the same origin and therefore the same slopes 
(see discussion in Sec. 4. 2). We then subtracted the fitted 
lines and bands from the original data-cube, pixel by pixel 
producing 'pure' continuum spectra. We visually inspected 
the subtracted data cube to ensure that no 'absorption' 
features due to overestimated lines/features were produced. 
From the 'pure' continuum spectrum, we built images 
in the ISOCAM LW5 (6.5-7.0 fim) and LW9 (14-16^m) 
filters. They represent the pure continua shortward and 
longward 10 /im, respectively. The NeII(12.7 /im) line was 
fitted with a Gaussian on the whole cube, even though we 
know that it is blended with the aromatic feature at 12.7 
/im which should be fitted with a Lorentzian. Therefore the 
Nell image contains also the 12.7 /im AFE. The resulting 
images are shown in Figures [3] and [H All emission shown 
hereafter are clipped at 3ct. 

Figure [jfl shows the 7.7 /im emission obtained through 
Lorentzian fits. Contours represent the SIV line emission 
obtained applying Gaussian fits (gray) and the 'pure' 15 
/im (LW9) continuum emission (black). In this image it is 
clear that the AFE is coming from a shell surrounding the 
HII region core and that this shell is particular bright in 
the PDR (i.e. in the North-East front). There are 2 bright 
compact knots on this shell (numbered in Figure [3]), not 
resolved at the 7.7/im resolution (^^7"). Most of the ionized 
gas emission arises from the cavity formed by the dust 
shell and it seems centered on a point like source. A similar 
result has been published recently by Zavagno et al. (|2006p 
on the massive galactic star forming region RCW 79. The 
15 /im pure continuum emission peaks close to knot#2 in 
the PDR dust shell and has a weaker emission in the 7.7 
fiui cavity where SIV peaks. Figure S] shows the 7.7 /tm 
image with the Nelll contours. The general picture is very 
similar to that shown in Figure [31 the ionized gas traced by 
the Nelll emission is mostly contained in the dust cavity 
and it does not show a secondary peak corresponding to 
the secondary peak seen in the SIV line emission, although 

^ Maps obtained from the CVF data have not been rotated 
north-south because at these wavelengths the images are un- 
dersampled and a rotation would conserve neither the flux nor 
the spatial distribution. To make the comparison among data at 
other wavelengths easier, we have rotated the Ha and the NIR 
images to the same ISOCAM orientation. 
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Fig. 3. Image: map from Lorentzian fit of the 7.7 nm dust 
emission feature. Gray and black contours: from map of 
Gaussian fit in the SIV line and 'pure' continuum in the 
LW9 filter (15/tm) emission respectively. 



this can be due in part to the poorer resolution at 15.6 
/tm than at 10.5 /tm. The main difference between the 
two ionized gas lines emission is that the Nelll emission 
is much more extended than the SIV emission towards 
the south part of the nebula. This is also the side where 
the HII region is partially bounded (Heydari-Malayeri and 
Lecavelier des Etangs lT994p and where the gas can expand. 
All other images in the other dust bands (not shown here) 
are very similar to the 7.7/(m emission. The NeII(-|-12.7 
/tm) distribution is also very similar to the distribution 
of the other AFEs everywhere but in the HII region core 
where there is a peak. This suggests, as expected, an 
increasing contribution of the Nell line to the 12.7 //m 
feature when approaching the exciting stars. We will come 
back to the relative contribution of the Nell line and the 
12.7 /im AFE in Sec. 4.1. 



3.3. J,H and Kg photometry. 

The J — H — Kg composite image of N4A (Figure [S|) shows 
an extended and diffuse emission whose morphology is very 
similar to the shape of the Ha emission seen in Figure[T]and 
with NIR color quite uniform across the nebula. Seven stars 
are visible (numbered in Figure [5l). Stars #1 and #3 are 
the bluest while the others are much redder. Star #4, at 
the north-west of the nebula, is bright and red. Star #2 is 
located at the border north of the nebula and shows dif- 
fuse emission around it. Photometry of all the seven stars 
is listed in Table 1 ordered by J magnitudes . 
The comparison with the position of the stars indicated by 
Heydari-Malayeri and Lecavelier des Etangs (|1994p as the 
two main exciting stars (marked as A and B in Figure [T]) 
suggests that they correspond to stars #1 and #3 in Figure 
[51 This association however, is not certain because of the 
lack of astrometry in their data (their Fig. 2). We believe 
that these are the exciting stars because if one associates 
the Ha knot visible in their Figure at about (0", -27") from 
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Fig. 4. Image: map from Lorentzian fit of the 7.7 //m dust 
emission feature. Contours: from map of Gaussian fit in the 
Nelll fine. 



the center of their field with our star#7, the position of the 
two exciting stars as indicated in their paper corresponds 
to the position of stars #1 and #3. 

The comparison between the band and the 15 /ini 'pure' 
continuum emission (shown in Figure [6|) clearly indicates 
that star#4 corresponds to the strongest peak of the pure 
continuum emission at 15 ^m. Figure [7] shows the Kg band 
emission contours on the SIV line image. The SIV emis- 
sion is almost entirely arising in the region where the 2 
bluest stars are located. Finally, Figure [8] shows that one 
of the AFE peaks (knot #2 in Figure [3|) is close but not 
totally corresponding to the bright continuum source (star 
#4) and that knot #1 is located north of the blue stars, 
very close to where is the red elongated feature visible in 
the NIR composite map (associated with #2 in Figure [5]). 
This suggests that this source and star #4 are embedded 
in a high density dusty region. 

Note however, that the uncertainty on ISO-CVF data as- 
trometry (<6") plus the fact that the images at differ- 
ent wavelengths are slightly displaced with respect to each 
other (by up 1 pixel over the whole ISOCAM spectral 
range) prevents us to use the ISOCAM astrometry to su- 
perpose images at different wavelengths. The method we 
used to align the CVF images with respect to each other 
and the CVF images on to the Kg band image, is matching 
the peak emission of the W01153 star (marked in Figure 
1)- 

From the comparison between the ISOCAM and NIR data 
we conclude that SIV and Nelll emission, i.e. the ionized 
gas, are mostly associated with the bluest stars in the cav- 
ity of the dust shell. The 15 /im continuum peak is mostly 
associated with the red star#4 very close but not coinci- 
dent to one of the peaks of the 7.7 emission (knot #2). 
Object #2 corresponds to knot #1 in the 7.7 /J,m emission. 



star IN b. 


T 

tJ mag 


(J-H) 


111 XI \ 


btar 1(A) 


13.63 


0.12 


0.24 


Star 2 


14.16 


0.32 


0.97 


Star 3(B) 


14.34 


0.07 


0.12 


Star 4 


14.50 


1.01 


1.26 


Star 5 


16.25 


0.93 


0.67 


Star 6 


16.26 


1.01 


0.71 


Star 7 


16.35 


1.04 


1.60 



Table 1. Observed colors of the 7 stars stars visible in NIR 
images. Numbers refer to Figure 5. Stars #1 and #3 are 
associated with stars A and B in Figure 1. 




Fig. 5. NIR Composite map oriented as the ISOCAM im- 
ages: red is Kg, green is H and blue is J. Point sources 
identified as stars are numbered as in Table 1. 



4. Discussion 

4.1. The ionized gas 

Figure [S] shows the Nelll/Nell ratio contours over the 7.7 
/im map. As explained above, this ratio has to be consid- 
ered a lower limit, because AFE at 12.7 /im and the Nell 
at 12.8 /tm are blended at the ISO spectral resolution. In 
general, the infrared Nelll/Nell ratio, practically not af- 
fected by extinction (Thornley et al. 120001 and references 
therein), is sensitive to the hardness of the UV radiation 
field. Thornley et al. "2000 studied how this ratio changes 
in a sample of star-forming regions belonging to starburst 
galaxies. They found that, excluding the two systems ( IZw 
40 and NGC5253) in their sample with low metallicity, the 
Nelll/Nell ranges from 0.05 to 1. In a comparable range of 
Nelll/Nell ratios, Brandl et al. (2006) have shown that in a 
sample of starburst galaxies the Equivalent Width (E.W.) 
of the 7.7 /tm feature does not vary significantly. At higher 
(> 1) Nelll/Nell ratios, a weak anti-correlation has been 
found by Wu et al. p006p in a sample of Blue Compact 
Dwarf galaxies, i.e. low metallicity systems. 

In N4 the Nelll/Nell ratio peaks in the Hll region 
core, right on top of the bluest stars, with a value equal 
to ~10 which is comparable to what has been found in IZw 
40 and in 30Dor (Thornley et al. 20Q0) and in other low 
metalhcity galaxies (Wu et al. [20061 Madden et al. 1 20061) . 
However, the global ratio is > 0.9, comparable with the 
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Fig. 6. Image: LW9 (15 fim) pure continuum emission; 
Contours= K band. 




Image= PAH at 7.7 jO. m 
Contours= K band emission 



Fig. 8. Image: PAH, 7.7 fim feature from Lorcntzian fit; 
Contours= K band. 




Image= SIV (10.5 jO. m) 
Contours— K band emission 



Fig. 7. Image: SIV at 10.5 /im from Gaussian fit; 
Contours^: K band. Contours are shown in two colors (black 
and white) for visualization purpose. 



typical value of the normal metallicity extragalactic HII re- 
gions. It is also noteworthy that there is an extended emis- 
sion in NeIII/NeII(+12.7 /im) ratio map in the south part 
of the HII region. 

Table 2 shows the total flux of the fine structure lines 
SIV, NeII(-hl2.7 /im) and Nelll. We checked whether the 
emission in the ionized gas lines of N4A agrees with the 
prediction of HII region models. We did not perform de- 
tailed calculation with models such as CLOUDY (Ferland 
et al. [T998J or similar, because we have only 3 lines avail- 
able, one of which is an upper limit. We instead compared 
our results with the model prediction from Stasihska p982p 
which gives the intensity of these lines relative to I(H/3). 
Heydari-Malayeri and Lecavelier des Etangs p994| give an 
1(11(3) intensity equal to 4.5x10"^^ erg/s/cm^ for a region 
centered on the exciting stars, of a radius equal to 21". In 
the same region we calculate the following ratios: 



IjSIV) 

I{NeIII) 
I{H(3) 

I{NeII) 
I{H(3) 



= 0.34 



0.6 



< 0.34 



(1) 



Assuming for N4 the physical parameters given in 
Heydari-Malayeri and Lecavelier des Etangs p994l their 
tables 6 and 8), i.e. Te//-40000-45000 K, n-100 
cm-3, [N]/[H ]=11.3e-6, [0]/[H]=2.6e-4, [Ne]/[H]=3.2e-5, 
[S]/[H]=6e-6, and a metallicity equal to Z/Z^oj^CS, we 
found that the corresponding model from Stasihska l|1982p 
(IDD2) predicts the following ratios: 



I{SIV) 

I{NeIII) 
I{H(3) 

I{NeII) 
I{HI3) 



= 0.4 



(2) 



0.6 



= 0.05 



These values compare very well with those observed in 
N4, listed in Eq. 1. The model also predicts an HII region 
radius of 4.2 pc, which corresponds to an angular radius 
of 16.5" (assuming a distance for LMC equal to 52 Kpc). 
This is quite the size of the Nelll/Nell emission if one does 
not include the extended emission on the south, outside the 
dust shell. Therefore, what we have called the HII region 
core, finds here its justification. It also roughly corresponds 
to the size of aromatic features cavity. 

We now compare the flux we obtained with previously 
published similar data from ISO-SWS. Giveon et al. (|2002p 
studied the mid-infrared fine structure lines emission of 
a sample of galactic and extragalactic HII regions, based 
on ISO-SWS spectra. N4A belongs to this sample. In or- 
der to compare our flux with those given in Giveon et al. 
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Aperture 


Fl0.5Mm(SIV) 


Fl2.8Mm(NeII+12.7) 


Fl5.6Mm(NeIII) 




IQ-is W/cm^ 


10~^^ W/cm^ 


10-18 W/cm^ 


Total 


68.7 


51.6 


65.0 


This paper in SWS aperture 


8.3 


6.8 


10.4 


Giveon et al. 2002 


5.7 


5.0 


13.7 



Table 2. Flux in the fine structure lines from ionized gas: first row gives the total fluxes; second row gives fluxes recovered 
in an aperture comparable to that of ISO-SWS. 



(|2002p . we should extract the flux in the exact same aper- 
ture as ISO-SWS (14"x20"). This is not possible with the 
ISOCAM data due to the fact that the pixel size of the 
images is 6" per pixel and any rebinning is unwarranted 
for quantitative calculation because the pixel size under- 
samples the PSF at all ISOCAM wavelengths. Therefore 
we can extract the flux in a box which is as close as pos- 
sible (12"xl8") to the size of the ISO-SWS aperture cen- 
tered on the HII region peak. The fluxes we recovered are 
listed in Table 2. They agree within 30% with the ISO-SWS 
fluxes. This is less than the total uncertainties, calculated 
taking into account the ISOCAM photometry uncertainties 
(~10%), uncertainties introduced by the fitting which we 
estimate to be ~10% and the ISO-SWS photometric un- 
certainties (~30%). Thus the total uncertainty is > 33% 
because we did not include in the calculation the uncer- 
tainties due to the differences in apertures and in pointings 
which cannot be estimated. 

The agreement between the fluxes we recovered from 
the ISOCAM data treated with the continuum plus feature 
fitting technique and the ISO-SWS fluxes, makes us 
confident on the reliability of the analysis presented here. 
An interesting point is that the ISOCAM NeII-hl2.7AFE 
flux accounts for all the flux in the ISO-SWS spectrum 
for Nell alone. This suggests that the contribution of the 
12.7 /im dust emission in the center of the HII region is 
negligible as it could have been guessed from the weakness 
of the other dust features. Therefore, in the HII region 
core, the Nelll/Nell value we calculated can be considered 
not an upper limit but the actual value which can now 
be directly compared with the typical ratios of poor 
metallicity systems given by Thornley et al. (|2000p . 



4.2. Dust emission 

4.2.1. Comparison with the current dust model 

The AFE carriers are thought to be Polycyclic Aromatic 
Hydrocarbons (PAHs) , composed by 50 up to few hundreds 
atoms, and transiently heated mainly by Far UV photons 
(6-13.6 eV). Their emission intensity is proportional to the 
product of the dust density and the ISRF, although it is 
certainly also depending on the hardness of the radiation 
field but in an unknown way. VSGs emit through a mecha- 
nism which is intermediate between the stochastic heating 
typical of PAHs and the thermal equilibrium of the classical 
big grains. In quiescent regions dominated by the general 
interstellar field, the dominant VSGs exciting mechanism 
is the stochastic heating which results in a smooth contin- 
uum peaking between 30-40 ^m. The emission mechanism 
becomes closer and closer to a standard thermal behavior 
when the radiation field gets more intense, such as in star- 
forming regions and PDRs. In these conditions, because of 
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Image= PAH at 7.7 ^ m 

Black Contoi.u-s= Neni/(NeII+12,7) ratio 




-30 O +30 
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(XINTEIl: R.A. 04 53 3.04 D£C -SS 56 EQUINOJO J^OOO 

Fig. 9. Image: PAH, 7.7 /im feature from Lorentian fit; 
Contours— NeIII/NeII(+12.7) ratio. Contour values are: 
0.7, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10.5. 



their quasi thermal behavior, the VSGs emission spectrum 
shifts toward shorter wavelengths. 

Our observations are in total agreement with the picture 
generally accepted to explain the dust feature and contin- 
uum emission at MIR wavelengths. First, we have shown 
in Sec. 3.2 that the morphology of the emissions in all dust 
features, including the NeII-|-12.7 /im if one exlude the HII 
region core, are very similar to one another, which strongly 
points towards a common origin and excitation mechanism 
for aU AFEs. 

Second, the spectra of N4 A presented in Figure[2]can be 
fully explained with the classical picture illustrated above: 
in the HII region core the continuum longward 10 /.tm in- 
creases just because VSGs reach the thermal equilibrium 
regime, get hotter, thus producing a shift of their thermal 
emission spectrum to wavelengths shorter than 30 ^m. The 
opposite happens in quiescent regions, such as that corre- 
sponding to the eastern CO peak outside the HII region, 
where in fact almost no continuum is detected. 

4.2.2. Origin of the continuum at MIR wavelengths 

Since we have built 'pure' continuum images at 6.75 and 
15 /xm, we can also investigate how these continua re- 
late to AFEs. The origin of the MIR continuum in star- 
forming regions, PDRs and diffuse ISM, is a controver- 
sial issue. In principle it can be due to hot VSGs which 
emission, as we have explained above, shifts into the MIR 
ISOCAM wavelengths and/or a continuum related to the 
AFE (maybe simply due to a superposition of unresolved 
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Fig. 10. Image: map from the 15 /im 'pure' continuum 
emission. Contours: 'pure' continuum at 6.7 fiia. smoothed 
to the 15 iitlh resolution. 



and weak bands). Figure [TOl shows the 6.7 /im 'pure' con- 
tinuum' contours on the 15 /im 'pure' continuum emission 
image smoothed at the same resolution. Maps have been 
aligned on the WOH53 star. The two continua are very 
similar suggesting a common origin. On the other hand 
their morphology is not totally similar to the dust feature 
emission morphology: they are more centrally peaked than 
the dust feature emission which forms a cavity. We there- 
for conclude that, at least in conditions similar to those of 
N4, i.e. in HII region complexes, most of the continuum at 
MIR wavelengths , even at wavelengths shorter than ''^10 
/xm, arises from VSGs and not from unresolved weak aro- 
matic bands. This might not be the case in more quiescent 
regions like diffuse ISM, as indicated by the spectrum cor- 
responding to the eastern CO peak (Figure[2]) where almost 
no continuum at all has been detected. 

4.2.3. Why is there no AFE in the HII region core? 

In this section we investigate the reason why no AFEs 
are detected in the HII region core. The depression of 
AFE in high ISRF environments has been already observed 
in the Milky Way (Cesarsky et al. , I1996bll99"6a|) . Small 
Magellanic Cloud (Contursi et al. , I2000p and in external 
galaxies (Galliano et al. 120031 Lu et al. 120031 Houck et al. 
I2004p . Recently, a high resolution study of the 3.3 /im AFE 
(not included in the ISOCAM-CVF wavelength range) with 
ISAAC at the VLT in NGC253 and NGC1808 (Tacconi- 
German et al. I2005p . has shown that AFE are strongly de- 
pressed with respect to the continuum close to where the 
super star clusters are located. For this reason these authors 
concluded that the PAHs may be a better tracer of B stars 
rather than more massive and earlier spectral type stars, as 
was already previously suggested by Peeters et al. (|2004p 
on the basis of the analysis of a wide sample of galactic HII 
regions observed with ISO. 

In the case of N4, almost no AFE is detected at all in 
the HII region core. This can in principle be explained with 
three possibilities: 



— 1) stellar winds are capable to evacuate the HII region 
from AFEs carriers but not from gas; 

— 2) Radiation pressure is capable to push away the dust 
but not the gas; 

— 3) AFEs carriers are destroyed in very extreme 
Radiation Fields (RF). 

Figure [3] shows that the 'pure' continuum emission at 15 fim 
has a secondary peak in the region where the ionized gas 
peaks. This is also visible in the MIR spectrum centered 
on SIV and Nelll peak (Figure [1] panel b) where almost 
no AFEs are visible but a significant continuum is still de- 
tectable. This shows that dust can survive in HII cores. 
As a consequence, the first two possibilities would imply 
that AFE carriers, VSGs and gas are decoupled which is 
physically quite improbable. We therefor conclude that the 
dominant mechanism has to be the destruction of AFE car- 
riers. 

The question now is whether dust destruction is due to the 
intensity or the hardness of the radiation field, or both. 
Recent results from Wu et al. (|2006p and Brandl et al. 
(|2006p have shown that PAHs E.W. stay quite constant for 
a wide range of UV hardness as traced by the Nelll/Nell 
ratio and start to slightly decrease only at high (> 1) 
Nelll/Nell ratio in low metallicity galaxies. More impor- 
tant is the result from Wu et al. (j2006p who showed that 
PAH E.W. is strongly anti correlated with the product of 
the radiation field intensity and hardness. 

We tried to investigate this issue as illustrated in Figure 
[TT] Here both panels show the 7.7 /xm dust feature emission, 
the 6.75 and 15 /im 'pure' continua emission as function of 
the ISRF energy density at 1600 A (left panel) and as func- 
tion of the Nelll/Nell ratio (right panel). The energy den- 
sity at 1600 A is in units of the solar neighborhood at the 
same wavelengths as measured from Gondhalekar Phillips 
and Wilson (1980) and it has been estimated in Contursi et 
al. 11998} ■ These authors calculated at each wavelength the 
average emissions in annuli centered on the main ionizing 
stars till 3 a level. The Nelll/Nell ratio is a tracer of the 
UV hardness. However, as already discussed in Sec. 4.1, at 
the ISOCAM spectral resolution, the Nell line emission at 
12.8 /xm is blended with the AFE at 12.7 ^m. Therefor, in 
the case of N4, this ratio is a lower limit everywhere but in 
the HII region core, where most of the flux can be ascribed 
to the Nell emission line alone, as it has been shown in Sec. 
4.1. In order to enhance differences among AFE and 'pure' 
continua emission, we have normalized each emissions to 
their maxima. 

Form Figure [Til we observe the following behaviors: 

— All dust components (AFEs, and the continua at the 
two wavelengths) show the same behavior: their emis- 
sion increase till a certain value going from the outside 
toward the HII region center. They reach a maximum 
presumably on the PDR and then sharply decrease in 
the HII region core. 

— The AFEs's drop is the steepest, followed by the the 6.7 
fim continuum and then by the 15 fiia continuum. 

— The AFE's drop occurs first, i.e. at ISRF density and 
UV hardness lower than for the pure continua emission 

— The decline of all emission appears steeper when plot- 
ted versus the Nelll/Nell ratio than when plotted as 
function of the ISRF energy density 
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These results definitely point to a scenario where, given 
the same conditions, AFE carriers are more easily destroyed 
than hot dust. In other words, these carriers are destroyed 
by UV photons less energetic than those necessary to de- 
stroy the carriers responsible for the 'pure' continuum emis- 
sion. If, following current dust models, these two types of 
carriers are associated with PAHs and VSGs, we can in- 
terpret this behavior as an evidence that PAHs and VSGs 
are both destroyed but with different efficiencies. Moreover, 
assuming that the VSGs population has a certain size distri- 
bution, we can reasonably think that the continuum emis- 
sion at lower wavelengths is dominated by smaller parti- 
cles than those responsible for the emission at longer wave- 
lengths. In this framework, the fact that the drop of the 
continuum at 6.7 /im is stepper than the drop of the 15 /im 
continuum, can be interpreted as the fact that the destruc- 
tion mechanism is more efficient on smaller molecule/grain. 
This means that: the destruction mechanism is likely to af- 
fect the dust size distribution. 

This conclusion may have some important implications for 
understanding the ISM physic in dwarf galaxies. Thornley 
et al. ((2000) found a Nelll/Nell ratio for IZw40 similar 
to what we measured in the very center of N4A (roughly 
a shell of ~10 pc in radius) but in a much larger region 
(~700x~1000 pc). This means that the phenomena of de- 
struction in such galaxies must occur in much larger re- 
gions than in normal galaxies. This can be explained by 
the fact that dwarf low metallicity galaxies have an aver- 
age UV spectrum harder than normal or starburst galaxies, 
due to a higher star formation rate per unit mass and to 
the presence of more young and massive star clusters. The 
UV photons can travel longer before being absorbed thanks 
to the fact that these systems have less dust (Madden et 
al. I2006p . All these factors point to a scenario where in 
dwarf galaxies, PAHs have a higher probability to be de- 
stroyed than in normal metallicity environments leading to 
an overall change of the galaxy dust size distribution. It 
seems also that, at least in not extremely poor metallicity 
systems, this phenomenon is more important than intrinsic 
chemical dust properties modification. 

Unfortunately, from this analysis no firm conclusion can be 
achieved about the possible different roles played by the 
ISRF intensity and hardness on the dust destruction mech- 
anism. However, a weak evidence that the latter has an im- 
pact stronger than the ISRF intensity might be deduced by 
the steeper drop all dust emission show with the Nelll/Nell 
ratio than that shown with the ISRF energy density. 

4.3. Comparison with hot dust properties in galactic 
environments. 

Metallicity is a parameter which may significantly affect 
dust. In low metallicity environments, such as those we 
expect to find in young distant galaxies, we expect dust to 
have intrinsic chemical-physical properties very different 
from those observed in normal metallicity galaxies, just 
because dust forms and grows in an ISM with very 
different chemical enrichment and mechanical inputs (from 
both stellar winds and supernovae). These might lead 
to un overall under-abundance of carbon based grains 
with respect to the abundance in normal metallicity 
environments but also to dust feature emission profiles, 
widths and intensity ratios different form those found in 
normal metallicity ISM. 



In Section 4.2.3 we have shown that in N4 dust destruction 
could be, at least in part, responsible for lack of small 
molecules/grains, already detected on much larger scale 
with IRAS in the Magellanic Clouds (Sauvage, Vigroux 
and Thuan [T990p . although we might have a total amount 
of carbon based molecule/grains lower than in Galactic 
HII regions. On the other hand, we have detected strong 
aromatic features very similar to what is observed in the 
Milky Way's ISM. In this section we will analyse in more 
details the N4 dust features widths and intensity ratios to 
understand their physical characteristics and how do they 
compare with similar observations in galactic environments. 



Dust feature intensities. 

Figure fT2l shows a comparison between the PDR spec- 
trum in N4 (panel c in Figure[2|) and the spectrum of one 
pixel in the PDR of NGC7023, after having subtracted 
a straight line continuum from both and scaled up the 
spectrum of NGC7023 to match by eye the 7.7 fj,m fea- 
ture of N4. There are no evidences, at least at the CVF 
spectral resolution (A/AA~50) of grain/molecules mod- 
ifications due to metallicity resulting in differences in 
the features profiles or peak wavelengths shifts. Indeed, 
the shape and the relative intensities of the 6.2, 7.7 and 
8.6 /im features in NGC7023 and N4 match almost per- 
fectly. There is an indication of a relatively lower 11.3 
fiui feature in N4 than in NGC7023, which may be due 
to a higher level of PAH ionisation in N4 . This is plau- 
sible, since N4 is more active than NGC7023, even in 
its PDR region as it is witnessed by the presence of ion- 
ized gas lines totally absent in the NGC7023 spectrum 
(but see more details below in this Section). This con- 
firms the result obtained by Vermeij et al. (2002) which 
found that the AFEs characteristics of LMC HII regions 
arc similar to those in the Milky Way. 
Dust features widths 

The analysis of the dust features conducted in this 
work is based on the assumption that their profiles are 
represented by Lorentzian. In this context, the widths 
of the features have a physical explanation. Boulanger 
et al. (|1998p . proposed that the emitting mechanisms 
of the AFE's carriers is the internal vibrational redis- 
tribution. The bands widths arise from intramolecular 
broadening processes in large molecules and are related 
to the very short life-time of the emitting levels. 
One of the byproducts of the fitting technique we 
have used, are the widths of the fitted Lorentzian 
through the entire data cube {i.e. wavelengths). We can 
therefore compare what obtained for N4 with the range 
of widths that Boulanger et al. fl998) found for typical 
galactic regions (namely p Ophiuchus and NGC7023). 
Figure fT3l shows such a comparison. We have calculated 
the mean widths of N4 obtained by the fit, taking into 
account all pixels with a signal greater than the 3 cr of 
the r.m.s. at each wavelengths (filled circles) and the 
median for the PDR region only, i.e. the same region 
shown in panel c of Figure [2] (filled triangles). The open 
stars corresponds to the mean FWHM values found 
for the two galactic regions by Boulanger et al. (jlQQSp . 
We exclude the 12.7 /im feature because we fitted it 
with a Gaussian and not with the Lorentzian as in 
Boulanger et al. 119981 The agreement is very good at 
all wavelengths showing that also this parameter is 
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Fig. 11. The 7.7 PAH emission (filled circle) 5/im and 15 /im pure continuum emission (triangles and diamonds 
respectively) versus the intensity of the radiation field at I6OOA (left panel) and Nelll/Nell ratio which traces the UV 
hardness (right panel). Values are above the 3a for the 3 wavelengths. In the right panel, the maximum at 7.7 /im and 
5 /im have been scaled to the maxima at 15 /xm to ease the relative comparison among the different emissions. 



very similar to that observed in galactic environments 
despite the fact that p Ophiuchus and NGC7023 are in 
ISRF much less intense than N4. 

Dust feature ratios 

Reach et al. (2000) found that the 11.3 /7.7 fj,m ratio 
was higher in SMCBl-1, a relatively quiescent cloud in 
the Small Magellanic Cloud (SMC) not associated with 
on-going star formation, than in various ISM regions in 
the Milky Way. This was attributed to an enhancement 
of the C-H bonds with respect to the C-C bonds, due 
to the higher H/C abundance ratio in the SMC than in 
the Milky Way. 

We have produced AFE feature ratio maps (not shown 
here), to study their variation across the region. The 
mean 11.3/7.7 /zm (C-H/C-C) value is 0.25 very 
similar to what found by Boulanger et. al. (1996) in 
p-Oph and smaller than the value in SMCBl-1 ( 0.8). 
A collection of 11.3/7.7 /im ratios in different galactic 
ISM by Lu et al. (1998) though, clearly indicates that 
this ratio has a large dispersion even in similar ISM 
phases. Therefore, the analysis of the dust feature ratio 
confirms that in N4 the dust feature ratio are compara- 
ble to what found in many other regions. An important 
result is that this ratio is quite constant on the dust 
shell, indicating that whatever is the state of the AFE 
carriers, i.e. ionized/neutral, hydrogenated or modified 
from metallicity, there are no strong differential effects 
on the dust feature emission, and therefore presum- 



ably on the dust conditions, at least at this spatial 
resolution (^ 2 pc assuming a LMC distance of 52 Kpc). 

The changes in the profiles of the dust features and 
their relative intensities can tell a lot about the physical 
state of AFEs carriers if associated with PAHs (Pauzat, 
Talbi and Ellinger Il997|, Le Page, Snow and Bierbaum 
[20031 Dartois an d D'H endecourt [19971 Bakes, Tielens 
and Bauschlicher I200T ) . While the first analysis is im- 
possible at the ISOCAM-CVF spectral resolution, we 
can analyze how AFE ratios change in N4 as function 
of the ISRF intensity. 

It is believed that PAHs in PDRs are mostly ionized. 
Ionization of PAHs leads to an enhancement of the in- 
tensities of the features arising from the C-C bonds (6.2 
and 7.7 /j,m bands in the ISOCAM wavelengths range) 
with respect to those arising from the C-H bonds (8.6 
/im in-plane and out-of plane (oop) solo, duo, trio and 
quatro at 11.3, 11.9, 12.7 and 14 /tm, respectively). 
Dehydrogenation weakens the C-H bonds. Therefore 
the ratio between any feature arising from C-H bonds 
with a feature arising from C-C bonds, should decrease 
as ionization and dehydrogenation increase. Figure [T4| 
upper panel shows the 11.3/7.7 /tni ratio (oop C-H/ vi- 
bration C-C) as function of the ISRF. There is a weak 
decreasing trend indicating that either ionization and or 
dehydrogenation are perhaps weakly increasing as ap- 
proaching the HII region core. On the other hand, in 
the very center there is an increase of this ratio proba- 
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A jdm 

Fig. 12. Comparison between the spectrum of the galactic region NGC7023 (thick soHd Hne) and the PDR spectrum of 
N4 (thin sohd hne), after continuum subtracted. NGC7023 spectrum has been scaled to match the 7,7 jim feature of N4. 



bly due to the fact that here smaller molecules (which 
emit at shorter wavelengths) are destroyed before larger 
species. Verstraete et al. (|1996p analyzed the oop C-H 
to in plane C-C dust features ratio, traced by the (11- 
13)/(6.2-|-7.7) /im ratio. Theoretical predictions show 
that this ratio is 2.5 for neutral PAHs and 0.25 for 
cations. Figure [T3] bottom panel shows a similar ratio 
where in place of the 11-13 emission features we calcu- 
lated the 11.3+12.8 yum emission. These ratios are well 
below those predicted for neutral PAHs everywhere in 
N4, and very close to the value for ionized PAHs, sug- 
gesting that PAHs are ionized in the whole region. Note 
that in the HII region core, this ratio is an upper limit 
because the 12.7 fim emission is almost entirely due to 
Nell rather than the AFE at 12.7 /im. We therefore con- 
clude that if PAHs are the carriers responsible for the 
dust features, in N4 they are almost totally ionized, and 
that their state of ionization and /or dehydrogenation is 
quite constant all over the region. Probably, in order to 
detect the transition between neutral to ionized PAH, 



one has trace the PAH emission from diffuse ISM to 
active regions. 

In conclusion, we have clearly shown that three fun- 
damental parameters of AFE, namely, emission intensities, 
features typical widths and emission ratios are all very sim- 
ilar to what is typically found in normal metallicity ISM. 
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Fig. 13. Typical widths of the AFEs obtained by the 
Lorentzian fit for all N4 (filled circles), the PDR only (solid 
triangles) and the typical values found with the same fit- 
ting technique by Boulanger et al. (1998) in p Oph and 
NGC7023 (open stars). 



[T5] This Figure includes for comparison the IR data ob- 
tained for two other LMC star forming regions, 30 Doradus 
(Rubio et al. 1998) and Nil (Barba et al. 2003). In each fig- 
ure's panel is plotted the zero age main sequence (ZAMS) 
adopted from Hanson, Howard and Conti ()1997p and the ex- 
tinction tracks for the reddening law of Rieke and Lebofsky 
p85|l . 

Note that the extinction for star #4 is likely to be much 
higher than the average as this star falls in the dust shell 
and on the 15 /im continuum peak. Four N4 IR sources show 
colors which are consistent with reddened main sequence O 
stars, while the other three are IR objects. 

From Figure [TSl we can conclude the following: 
- 1) Star #1 and #3 are probably moderately reddened 
main sequence O stars with Av ~ 3. They are the main ion- 
izing stars responsible for the HII region (Heydari-Malayeri 
and Lecavelier des Etangs [l994p and labeled as A and B in 
Figure 1. 

-2) Stars ^5 and ^6 arc highly reddened main sequence 
O stars with Av > 10. They do not have optical counter- 
parts. 

-3) The IR objects labeled #4, #5, and #2 have large 
IR excess and their colors are not that of reddened young 
main sequence stars. They lie in the border of the molec- 
ular cloud seen in CO (Contursi et al. I1998P and at the 
edge of the cavity formed by the ionizing stars. These three 
objects have J band brightness corresponding to Class I 
and /or Herbig Ae/Be objects candidates (see Brandner et 
al. I200ip . Stars #7 and #4 have the largest IR excess of 
the sample. They have IR characteristics of massive young 
Stellar Objects(YSO's) similar to N122 in 30Dor. Object 
^2 is an extended object which shows nebulosity around 
and a second red component. 
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Fig. 14. Different dust feature ratios as function of the 
ISRF. Upper panel: 11.3/7.7 ^vci dust feature ratio, which 
traces the [out of plane C-H]/ [vibration C-C] ratio. Bottom 
panel: (ll-13)/(6. 2-1-7. 7)/xm dust feature ratio which traces 
the [out of plane C-H] /[in plane C-C] dust features ratio. 



5. The stellar content 

The NIR data presented in this paper unveil at least 7 
bright sources. Figure [T51 show a NIR color-magnitude di- 
agram (right panel) and JHKg color-color (left panel) of 
stars numbered as in Figure 5 whose photometry is reported 
in Table 1. These are shown as red pentagons in the Figure 



5.1. The nature of star ^ 4 

This is the brightest member of a small cluster and it is of 
particular interest because it corresponds to a bright point 
source emission in the hot dust pure continuum emission. 
Its position in the NIR color-color and color-magnitude di- 
agram is similar to that of IR sources found in 30 Doradus, 
i.e. source W30, N122 and N126 (Rubio et al. 1998). These 
30 Doradus sources were resolved into multiple compact 
young massive systems from high resolution HST/NICMOS 
observations, with a dominant IR source (Walborn, et al 
2002, Walborn, Maiz-Apellaniz, Barba, 1999). The IR color 
of star #4 can be either interpreted as early O type star 
with Av= 20 and with a IR K excess of about 4 magnitudes 
or a multiple object containing several ZAMS 03 stars ex- 
tinguished by 20 magnitudes. In fact, the star profile is 
wider than the point spread function derived for the stars 
in the Ks image, which could be a hint of the presence of a 
much redder component, thus favoring the hypothesis of a 
multiple system. 

We can check the results on the nature of star #4 
obtained from NIR photometry, with the comparison of 
our ISOCAM and NIR data with those published by the 
ISOGAL team (Felh et al. [2002]) . The LW2 and LW3 
broad band magnitudes obtained through PSF fitting are: 
m(LW2)= 6.1 and m(LW3)= 3.4. Then, the position of star 
#4 on the [15]-[7]-[15] magnitude-color diagram published 
by Felli et al. (|2C)02p falls exactly in the region occupied by 
Young Stellar Objects (YSOs). 
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Fig. 15. NIR color-magnitude (left) and color-color (right) diagrams of N4 with the IR sources showed as red filled 
points. For comparison we include the 30 Doradus IR sources as open circles and the Nil IR sources as open squares. 
In the color-magnitude diagram the upper zero-age main sequence between 03 V and 09 V corresponds to a distance 
modulus of 18.6 and is indicated with a dashed line. The reddening track for a normal 03 V star is plotted with a dotted 
line and extends to Ay = 20 mag. In the color-color diagram, the main-sequence locus from 03 V to M2 V and the 
cool-giant branch are indicated by a solid line. The reddening tracks for normal 03 V and cool giant stars are plotted as 
dashed lines, with crosses indicating Ay = 10 and 20 mag. Sources 1, 4 and 7 show IR excesses which cannot be explain 
from reddening. 



A third indication that this objects contains indeed an 
YSO, is given by its CVF spectrum shown in Figure [T6l 
It shows a deep silicate absorption at ~ 10 /im which sug- 
gests this is a high embedded source. This is also confirmed 
by the fact that star #4 corresponds to the strongest peak 
of the 15 /im pure continuum emission. In the spectrum 
are also visible two weak absorption features at 13 fim 
~15/Ltm. We checked whether they are due to memory ef- 
fects of the detector after a glitch. The result is that most 
likely, the feature at 13 /im is fake but that at 15 /im seems 
real. If it is so, it could be associated with the solid CO2 
bending mode at 15.2 /xm or to gas phase CO2 at 15.0 fim. 
Both types of emission are detected in YSOs (van Dishoeck 
[2004). The fact that the overall spectrum shown in Figure 
\W\ docs not look like a typical YSO MIR spectrum (van 
Dishoeck I2004p could be due to the ISO spatial resolu- 
tion; the ISOCAM beam is in fact sufficiently large (2 pc 
at the LMC distance) to include a considerable amount of 
ISM surrounding the YSO, as suggested by the presence of 
AFEs. To our knowledge this would be the first time that 
the CO2 absorption band is (marginally) detected in extra- 
galactic YSOs, although SPITZER will certainly discover 
and study many more. Since such feature, together with the 
solid CO2 stretch mode at 4.25 /im and II2O at 3 fim (out 
of our wavelength range) has been observed at higher res- 
olution with ISO-SWS in some Ultra Compact HII regions 
(Roelfsema et al. 11998 !) this would support our hypothesis 
that star#4 is an embedded YSO. 
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Fig. 16. ISOCAM CVF spectrum of star #4 corresponding 
to the 15 /im pure continuum peak. 
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6. Conclusion 

In this paper wo have analyzed the ISOCAM CVF and new 
JHKs photometry data of the HII region complex N4 in the 
Large Magellanic Clouds. The analysis has two principal 
aims: 1) to relate dust, gas and stellar content; 2) to look 
for effects of the lower than galactic metallicity on the dust 
features. The main conclusions are the following. 

— The overall MIR spectral characteristics of N4 can be 
fully explained with the current dust models. The con- 
tinuum + dust features + gas lines properties are those 
of a typical galactic HII complex which can be inter- 
preted as the combination of three primary components: 
HII region core+ PDR+ parental molecular cloud. 

— Thanks to a fitting technique through the entire CVF 
data cube, we can produce maps in the single dust fea- 
tures, gas lines and pure continuum. We find that the 
HII region core is completely devoid of dust features. 
This dust cavity contains most of the ionized gas. The 
pure continuum emission seems to have the same ori- 
gin at all wavelengths, but it is probably arising from 
a population of grains different from the AFE carries. 
Following the current dust models this component can 
be totally ascribed to VSGs. 

— The comparison between the ISM emission and the stel- 
lar content of N4 derived from NIR photometric data, 
reveals that the two bluest stars are the main ioniz- 
ing stars of N4A: they reside in the HII region core, 
i.e. in the center of the AFE cavity and where the ion- 
ized gas peaks. The pure continuum emission peaks on 
a point source which shows strong evidences to be a 
group of massive young stars containing a deeply em- 
bedded YSOs. 

— A detailed comparison of the relative spatial distribu- 
tion and intensities of the dust features as function of 

the intensity of the radiation field, suggests that the HII 
region core does not contain PAHs because they are 
destroyed. We show that this mechanisms is selective 
depending on the grains sizes, the smallest grains/big- 
molecules, being destroyed first. This can lead to an 
overall modification of the grains size distribution which 
can be significant in relatively low metallicity dwarf 
galaxies, where these processes act on large portions of 
their ISM. We also argue that this mechanism is more 
important than any intrinsic dust modification due to 
the relatively low metallicity of LMC. 

— The comparison between the widths and the intensity 
ratios of the dust features in N4 and in other galactic 
environments shows that there are no evident and sig- 
nificant signs of metallicity effects on AFE's carriers in 
N4 

— The analysis of the dust features intensities ratios and 
their comparison with model predictions, indicates that, 
if PAHs are the particles responsible for the AFEs, they 
are fully ionized over the entire analyzed region. 

— The IR stellar content of the N4 shows several bright 
IR sources with characteristics of massive early O type 
stars similar to those found in SODoradus. IR spec- 
troscopy of these sources would confirm their very young 
and massive nature. 
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